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Abstract 

A two-dimensional cellular automaton model of traffic flow with open bound- 
aries are investigated by computer simulations. The outflow of cars from the 
system and the average velocity are investigated. The time sequences of the 
outflow and average velocity have flicker noises in a jamming phase. The low 

density behavior are discussed with simple jam-free approximation. 
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I. INTRODUCTION 



Models of traffic flow have relations to wide varieties of physical systems. A traffic flow 
system is one of asymmetric exclusion processes. They are non-trivial statistical mechanical 
systems because of lack of detailed balance. Studies of these non-trivial system explore 
the profound structure of statistical mechanics. Studies of exclusively interacting particle 
systems like traffic flow also relates to equilibrium and non-equilibrium properties of granular 
flows, surface growth, dynamics of defects in solids, and so on. The model I will discuss here 
may be one of the simplest examples of non-equilibrium colliding granular flows. 

Traffic flow problems have been studied mainly through fluid dynamics, car-following 
models, coupled map lattice models, and cellular automaton (CA) models. Many attempts 
have been made to apply CA modeling to complex phenomena including fluid because of 
computational simplicity. Cellular automaton modeling of traffic flow is one of the recently 
developing area. One of the simplest CA models of traffic flow in one-way expressway is 
the rule-184 elementary CA 0, which is a simple asymmetric exclusion rule. In spite of 
the simplicity of the model, it shows a phase transition from a freely moving phase at low 
vehicle density to a jamming phase at high vehicle density. The computational simplicity of 
CA models also enables us to take many realistic features of traffic problems into account. 
More realistic models considering speed variation of cars or effects of blockades have been 
investigated in one-dimensional models 1/f fluctuation has been observed in both 

actual expressways || and models |||7|]. Self organized criticality has also been studied ||. 

Traffic networks, for example a traffic system of a whole city or an expressway network, 
consist of many complicated ingredients. It is very hard to model the whole features of 
traffic networks. Two-dimensional CA models of traffic flow, therefore, are very abstract 
models of traffic networks. One of the simplest two-dimensional CA model of traffic flow 
has been investigated by Biham, Middleton and Levine (BML) ||. Their model is a simple 
extension of the rule-184 CA to two dimension. Cars are distributed on a square lattice of 
N x N sites with periodic boundaries with both in the horizontal and vertical directions. 
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They found a sharp transition between a freely moving phase at low vehicle density and a 
jamming phase at high vehicle density The characteristics of the transition were studied 
by Nagatani [I0j and Fukui and Ishibashi flIT] , |T2|j . Two type of jam phase was discussed 



by Tadaki and Kikuchi [|13|,|T4[]. The model has been extended to take into account the 
probability of changes in vehicle directions ]T5| , P^ ]. 

Two-dimensional CA models of traffic flow show many physically interesting phenomena, 
phase transitions and self organization. Cellular automaton modeling of traffic systems, 
however, is a toy model. It should be clarified what features strongly depend on the model 
itself. Characteristic features of the BML model are, for example, deterministic dynamics, 
periodic boundaries, restrictions on the car destination, road arrangement without traffic 
queue [|17| and so on. In this paper, open boundary conditions instead of periodic ones are 
used to investigate the emergence of traffic jam in two-dimensional CA model. 

The organization of this paper is as follows: The model is given in 0. The dynamics 
is described with binary arrays. The outflow of cars from the system is investigated in [TTI] . 
The jam-free approximation is discussed. In [TV] the average velocity of cars is investigated. 
Section [V] is devoted for discussions. 



II. MODEL 

The model is the same as the model-I of BML except the boundary conditions. On the 
contrary to the original BML model, cars are injected probabilistically both the left and 
lower boundaries of the system and flow out deterministically from both the right and upper 
boundaries. 

Up-directed and right-directed cars are exclusively distributed in N x N square lattice. 
Each site is empty or occupied by one up-directed or right-directed car. Cars can move one 
step at a time if and only if the adjacent site in the destination is empty. There is a traffic 
light controlling the whole system as up-directed cars can move only at even time steps and 
right-directed cars can only at odd time steps. 
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The number of right-directed (up-directed) cars at a time t and a position r = 
(1 < h j < N) is expressed by a binary array /v(t) = {0, 1} (v?(t) = {0, 1}) [0. The bulk 
dynamics, namely dynamics for bulk sites (1 < i < N, 1 < j < N for up-directed and 
1 < z < iV, 1 < j < iV for right-directed cars) can be expressed as 

/J?(t + 1) = a(t)jlf(t) {fir+x(t) + W+s(t)} 

+a(t) {1 - fi P (t)} {1 - w{t)} lirs(t) 

+ {l-a(t)}Mt), (1) 
u^t + 1) = {1 - <r(t)} ^(t) {^ +s/ (t) + z/ r - +?7 (t)} 

+ {i - (7(0} {i - M*)} {i - ^(t)}^,_. ? (t) 

+ff(^fW, (2) 

where x and y denote unit vectors of right and up directions respectively. The binary function 
a(t) = t mod 2 represents the control by the traffic light. The condition n^(t)i>^(t) = holds 
because one site can not be occupied both with up and right cars simultaneously. 

The first term a(t)fi^(t) {/z^ + ^(t) + u^{t)} in eq. ([]]) denotes that a right-directed car 
remains at the site r if the right adjacent site is occupied by a right-directed or up-directed 
car. The injection of a right-directed car from the left adjacent site is given by the second 
term a(t) {1 — /v(£)} {1 — Vr(t)} Hr~x{t)- The last term {1 — o~(t)} /v(^) shows that a right- 
directed car does not move at odd time steps. The same decomposition of eq. @ can be 
done for the dynamics of up-directed cars. 

Cars are injected from lower and left sides of the system. If the site on the edges of the 
system is empty, a car is injected with a probability p. The injection of right-directed cars 
on the left edge r = (1, j) (1 < j < N) is given by replacing the injection term (second 
term) in eq. (|l|) with probabilistic injection: 

fif(t + 1) = a(t)/If(t) {fir+x(t) + W + x(t)} 

+a(t){l-Mt)}{l-Mt)}f(p) 

+ {l-a(t)}ii r it), (3) 
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where f(p) = {0, 1} is a function which returns unity with a probability p. The injection of 
up-directed cars on the lower edge r = (i, 1) (1 < i < N) is given by 

Vr{t + 1) = {1 - (j(t)} V?(t) {n? + y{t) + Vf+y{t)} 

+ {1 - a(t)} {1 - Mt)} {1 - ^(*)} f(p) 

+a{t)v?{t). (4) 

Cars flow out from the upper and right edges of the system deterministically. The 
dynamical equations for cars on the upper and rights edge are given by deleting the first 
terms in eqs. ([!]) and (§). For sites on the right edge r — (N,j) (1 < j < N) 

fj^t + 1) = a(t) {1 - firit)} {1 - v F {t)} nr_g(t) 

+ {1- a(t)} fif{t) (5) 

gives the dynamics of right-directed cars on the right edge. For sites on the upper edge 
f = (i,N) (l<i<N) 

V ? {t + 1) = {1 - a(t)} {1 - flfit)} {1 - Mt)} "r-y(t) 

+a{t)w{t) (6) 

gives the dynamics of up-directed cars on the upper edge. 

In the current simulations, the system has no car at the initial time t = 0. Cars are 
injected with eqs. (|3]) and probabilistically and run deterministically obeying eqs. ([[]) 
and (||). If cars reach the edges of the system, they flow out by eqs. (|5]) and (H). In the early 
traffic light cycles to O(N), the front lines of right-directed and up-directed cars collide to 
form a global traffic jam configuration in case p > p c (p c is discussed later). The global jam 
is sorted out with the maximum throughput, where the number of cars per site is p = 2/3. 
Then new small jam clusters are created and sorted out again and again. Figure [1] shows a 
snapshot of the system. In the simulation, the system runs 200 N times (lOOiV traffic light 
cycles) from the initial condition for relaxation. And quantities discussed later are observed 
for 200iV < t < 400iV. 



III. OUTFLOW 



The first quantity we observe is the outflow of cars from the system. By virtue of the 
dynamics, eqs. ([5|) and (|6|), the outflow is the number of cars appearing on the upper and 
right edges of the system. The average outflow p out is defined as the average number of cars 
appearing on the edges per site and traffic light cycle. The results of the simulation is given 
in Fig. |. 

In the low injection p region, the system can be assumed to be free from traffic jam. In 
this case the injection process will be controlled only by the number of cars which stay on 
the lower and left edges of the system. The number of up-directed (right-directed) cars in 
each column (row) is given by p on tN. There are 2p ont N cars on the left edge of the system. 
These cars prevent the car injection from the left side of the system, and the remaining 
(1 — 2p out )N sites can accept the car injection. At the next time step, therefore, the number 
of cars injected on the left edge will be (1 — 2p out )Np. The equilibrium condition between 
the injection and the outflow 

Pout = (1 - 2p ont )p (7) 

gives 

Pout = ^ (8) 

as the average outflow. This naive estimation of the outflow (jam- free approximation) agrees 
well with the simulation results for p < p c (p c ~ 0.2). The extrapolation of eq. (||) to p = 1 
gives p ou t = 1/3, which corresponds to the maximum throughput. 

In the high injection p region, the system has traffic jam clusters in the bulk area. The 
outflow pout is suppressed and lower than that given by the jam- free approximation discussed 
above. 

The time dependent behavior of the outflow p ut(^') (f denotes the traffic light cycle and 
t' — 0, . . . , T — 1) and its power spectrum 
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t'=0 



is observed, where T is the maximum traffic light cycles obeying T = 2 T < lOOiV [ll|. In 
the high p region (Fig. f|), the power spectrum of p ut(t) shows Ik ~ kr a behavior (Fig. ^|). 
This shows the existence of the self-organized jam clusters in the bulk system. In the low 
p region (Fig. |3|), on the other hand, p ut(t') shows random fluctuation around the average 
and the power spectrum is beared with week flicker noise. 



IV. AVERAGE VELOCITY 

The average velocity of the cars is the number of cars moving during one traffic light 
cycle (namely two time steps). The arrays \i and v are binary ones. Thus the average 
velocity, naively saying, is half of the Hamming distance between fi(t) and fi(t + 2) (v(t) 
and v(t + 2)). Caution must be payed to tread the edges of the system. For example 
for right-directed cars, if the left edge site is empty the site must be excluded from the 
calculation of the Hamming distance because of the probabilistic injection. If the right edge 
site is occupied by right-directed car, the site must be excluded because of the deterministic 
outflow. The same treatment is applied to up-directed cars. Figure || shows the results of 
the simulation. 

In the low p region, the jam-free approximation gives the outflow p ou t as discussed in PH . 
There are 2p ont N 2 cars in the system. By the assumption of no jam in the bulk area, cars 
are distributed randomly. There are pl nt N 2 colliding pairs of cars. The number of freely 
moving cars will be 2p out N 2 (l — (l/2)p out ). And the average velocity is 

v = l--—^—. (10) 
2 1 + 2p K J 

This estimation agrees with the results of the simulation less than the case of the outflow. 
The discrepancy seems to come from the effect of collisions with more than two cars. These 
effects are expected to decrease faster than two-car collision in large systems. The simulation 
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results seems to show that the jam-free approximation becomes better with increasing the 
system size. 

At the critical injection p c ~ 0.2, the average velocity shows a sharp phase transition 
with its sudden decrement by the formation of jam clusters. It increases gradually with the 
increment of p above p c . The behavior of v just above p c shows strong finite size effects. 

The time dependent average velocity v(t') and its power spectrum 



shows the same characteristics as those of the outflow (Figs. |7| and |8|). The high p case shows 
the flicker noise as J k ~ k" 13 with (3 ~ 1.2 (Fig. |) reflecting the emergence of traffic jam. 
The low p average velocity seems to have week flicker noises. This shows the temporary 
formations of small traffic jam. 



In this paper a two-dimensional cellular automaton traffic flow model with open bound- 
aries was investigated by computer simulation. The bulk dynamics is deterministic. Cars 
are probabilistically injected from the left and lower sides of the system and flow out from 
the right and upper sides deterministically. 

The average outflow pout, which is the number of cars flowing out from the right and 
upper sides per traffic light cycle and per site, obeys p out — p/(l + 2p) in the low injection 
region (p < p c ), where p is the injection rate. This is well understood with the jam-free 
approximation. High injection p > p c causes the emergence of traffic jam clusters in the 
system and suppress the outflow. In the high injection region, the time-dependent behavior 
of the outflow shows flicker noises. 

The average velocity of cars was also investigated. The jam-free approximation value 
v = 1 —p/(l + 2p) /2 does not well agree the simulation results. The reason seems to be the 
many car collision effect which will be suppressed in large system size. The average velocity 
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V. DISCUSSION 



s 



shows sharp phase transition at p c ~ 0.2. In the high injection region, the average velocity 
is suppressed by the emergence of traffic jam clusters. The time-dependent behavior of the 
average velocity also shows flicker noises in the high injection region. In the low injection 
region, it shows week flicker noise because of the formation of temporal jam clusters. 

In the original BML model, which has periodic boundaries, cars are freely moving with 
v = 1 in the low density region. On the contrary, the current system with open boundaries 
shows the average velocity less than unity even in the jam- free state. The reason of the 
difference is as follows:As well known in one-dimensional case (Wolfram's rule-184), temporal 
jam clusters in the low density region are sorted out and form the maximum throughput 
current with the local density p = 1/2. Sorting out of jam clusters gives the average velocity 
v = 1. In two-dimensional cases, temporal jam clusters are also sorted out and form the 
maximum throughput current with p = 2/3. In the periodic boundary cases, once the 
coherent maximum throughput currents are created, they dominated the whole system and 
new traffic jam clusters are hardly created. In the open boundary cases, on the other hand, 
the coherent maximum throughput currents flow out from the system and the incoherent 
currents are injected. These incoherent injections form new traffic jam clusters and suppress 
the average velocity. 

In this paper, I called the event at p ~ 0.2 as phase transition. The event is not a 
phase transition in the strict sense of the word. No critical behavior is found at the point. 
The value of the average velocity shows sharp discontinuity at p ~ 0.2. It has finite value 
above the point however. An adequate order parameter is needed to strictly define the phase 
transition. 

In Fig. ^| the average velocity increases above the point p ~ 0.4. As mentioned above, 
the maximum flow with p = 2/3 is formed behind jam clusters. The contribution of the 
maximum flow to the average velocity is expected to grow with the injection rate p. The 
increment of p also seems to enhance the deterministic feature of the injection process. 
These factors may contribute to the increment of the average velocity. On the contrary, the 
increment of the injection rate p contributes the formation of jam clusters which decreases 
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the average velocity. The competition of this cluster formation effect and the previous two 
factors decides the behavior of the average velocity. The investigations of the statistical and 
dynamical properties of spatial structures may clarify the behavior of the system. 

Tadaki and Kikuchi shows the existence of two types of jam phases in the BML model 
|TB|jri| ]. The current model with open boundaries seems to have only one jam phase. The 
high density random jam phase found in periodic boundary case seems to be one of the 
finite size effects. In viewpoints of statistical mechanics, finite size effects will be neglected 
in realistic macroscopic systems. On the contrary, real traffic network systems are finite and 
the finiteness may be important factor of the system. Observations of real traffic network 
systems are expected. 
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FIGURES 

FIG. 1. A snapshot of the simulation. The system size is 100 x 100 and p = 0.4. The black and 
gray dots show right-directed and up-directed cars respectively. There are some local jam clusters. 
Jam clusters are sorted out with the maximum throughput (p = 2/3) . 

FIG. 2. Average outflow from the system for the system size 50 x 50, 100 x 100, 200 x 200 and 
400 x 400. The bold line gives jam-free approximation p out = p/(l + 2p) which is given in the text. 

FIG. 3. Time dependent behavior of p out (t) for p = 0.1 with the system size 400 x 400 and its 
power spectrum 1^. The bold line is fitted with the power spectrum 10 1 < k < 10 3 by the method 
of least square. 

FIG. 4. Time dependent behavior of p out (t) for p = 0.6 with the system size 400 x 400 and its 
power spectrum 1^. The bold line is fitted with the power spectrum 10 1 < k < 10 3 by the method 
of least square. 

FIG. 5. The behavior of the exponent a. The exponents are calculated to fit the power spectrum 
Ifc within 10 1 < k < 10 3 with the method of lease square. In the high p > p c region, the exponent 
behaves constant a ~ 0.8 which depends on the system size. Below the critical p c the exponent 
sharply decrease because the low p power spectrum shows week flicker noise. 

FIG. 6. Average velocity for the system size 50 x 50, 100 x 100, 200 x 200 and 400 x 400. The 
bold line is v = 1 - (l/2)p/(l + 2p). 

FIG. 7. Time dependent behavior of v(t) for p = 0.1 with the system size 400 x 400 and its 
power spectrum J^. The bold line is fitted with the power spectrum 10 1 < k < 10 3 by the method 
of least square. 

FIG. 8. Time dependent behavior of v(t) for p = 0.6 with the system size 400 x 400 and its 
power spectrum J^. The bold line is fitted with the power spectrum 10 1 < k < 10 3 by the method 
of least square. 
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FIG. 9. The behavior of the exponent (3. The exponents are calculated to fit the power spectrum 
Jfc within 10 1 < k < 10 3 with the method of lease square. In the high p > p c region, the exponent 
behaves constant (3 ~ 1.2 which depends on the system size. Below the critical p c the exponent 
sharply decrease because the low p power spectrum shows week flicker noise. 
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